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On convergence and stability 
conditions of homotopy perturbation 
method for an inverse heat conduction 
problem 


Q. Jannati and A. Zakeri* 


Abstract 


In this paper, we investigate the application of the Homotopy Pertur- 
bation Method (HPM) for solving a one-dimensional nonlinear inverse heat 
conduction problem. In this problem the thermal conductivity term is a lin- 
ear function with respect to unknown heat temperature in bounded interval. 
Furthermore, the temperature histories are unknown at the end point of the 
interval. This problem is ill-posed. So, using the finite difference scheme and 
discretizing the time interval, the partial differential equation is reduced into 
a System of Nonlinear Ordinary Differential Equations (SNODE’s). Then, 
using HPM, the approximated solution of the obtained Ordinary Differen- 
tial Equation (ODE) system is determined. In the sequel, the stability and 
convergence conditions of the proposed method are investigated. Finally, an 
upper bound of the error is provided. 


Keywords: Homotopy perturbation method; Diffusion equation; Disceretiz- 
ing method; Inverse problem. 


1 Introduction 


Inverse heat conduction problems are used to describe many important phe- 
nomena in physics, chemistry, mechanics, etc. There has been a great amount 
of investigation to solving inverse heat conduction problems in one and multi 
dimensional spaces. Many effective methods have been provided. However, 
lots of inverse heat transfer problems, which arise in natural phenomena, 
such as radiational heat transfer, modelling of case hardening, gravimetry, 
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and etc., have nonlinear forms and so they are not solvable with analytical 
methods, but unfortunately, most of presented methods are useful just for 
solving linear forms. 

Usually these problems are ill-posed in the sense of Hadamard. There- 
fore, the regularization method is a successive technique for solving ill-posed 
problems and it may be applied to entire class of problems which arise from 
physical observations. 

Beck et al. have investigated an inverse problem in one-dimensional space 
with two general procedures, function specification and regularization meth- 
ods, and a method of combining these, trial function method, and have im- 
plemented all of these methods in a sequential manner [3]. Lesnic et al. 
in [15] and [19] have considered a special case of distributed (identification) 
parameter problems in one-dimensional spaces. They have shown that for 
a one-dimensional quasi-heterogeneous material with square-root harmonic 
conductivity, a single measurement of the conductivity and the flux on the 
boundary is sufficient to determine uniquely the unknown physical parame- 
ters and the solution function. Alivanov considered the solution of inverse 
problems by analytical approaches [2]. Qu and Dou [20], Lewandowski [18], 
and Jia et al. [16], have studied the nonlinear diffusion equation and provided 
some numerical techniques. Shidfar and Zakeri in [21] - [23] have investigated 
the existence and uniqueness of a solution for a two-dimensional nonlinear 
inverse diffusion problem. Also, Zakeri et al. have begun their research by 
a Cauchy inverse problem and found a solution by HMP method [26] and 
in continuation they have gone on by an inverse heat conduction problem 
and solved it by the HPM again [27]. Also, they have applied an approach 
which contained a difference method and the HPM together, and solved the 
problem with a reliable accuracy [28]. 

In continuation of above researches our intend is investigation of suffi- 
ciently condition for HPM for solving inverse heat conduction problems. 

In next section, the HPM is introduced shortly and in Section 3, an ap- 
proximated solution for the inverse heat conduction problem is obtained via 
HPM . Then, the stability and convergence of the above mentioned method 
are studied in Section 3 and Section 4, respectively . Some numerical results 
are illustrated by some tables and figures in Section 6. Finally conclusions 
and some suggestions for more research are given in last section. 


2 Basic concepts of HPM 


In this section we introduce the basic concepts of the HPM , according to [26], 
in brief. 
Consider the following nonlinear equation 


A(u) — f(r) =0, reEQ, (1) 
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with the boundary conditions 
B(u, 0u/On) =0, rel, 


where A is a general differential operator, B is a boundary operator, f(r) 
is a known analytic function and I is the boundary of the domain 2. The 
operator A can be generally divided into two parts L and N, where they are 
the linear and nonlinear parts of A, respectively. So, Eq. (1) converts into 
the following form 


L(u) + N(u) — f(r) = 0. (2) 
In [10], He constructed a homotopy H : © x [0,1] —+ R which satisfies: 
H(v, p) = (1 — p)[L(v) — L(vo)] + pP[A(v) — f(r)] = 9, (3) 
or 
H(v,p) = L(v) — L(vo) + pL(vo) + p[N(v) — f(r)] = 9, (4) 


where r € Q, and p € [0,1] is called the homotopy parameter, and vp, is an 
initial approximation for the solution of Eq. (1) which satisfies the boundary 
conditions. Consequently 


H(v,0) = L(v) — L(vo) = 9, 
H(v,1) = A(v) — f(r) =0. 


Now, when p varies from 0 to 1, the homotopy H(v,p), changes from 
L(v) — L(vo) to A(v) — f(r). 

Applying the perturbation technique due to the fact that 0 << p <1is 
considered as a small parameter, we can assume that the solution of Eq. (3) 
or Eq.(4) can be expressed as a series in the form 


v=vp + pu, + p2ve +... : 


When p —> 1, Eq.(3) or Eq.(4) corresponds to Eq.(2) and so v becomes the 
approximate solution of Eq. (2) i.e; 


u=limv=vup9+utug+.... (5) 
pol 


The series (7) is convergent for most cases and the rate of convergence de- 
pends on A(v) [11,12]. 

In the next section, a nonlinear inverse heat conduction problem is consid- 
ered. We discretize the time interval by means of backward finite difference 
method and apply the HPM. Then, the approximate solution of the problem 
is yield. 
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3 Solution of nonlinear inverse heat conduction problem 
by HPM 


Let T > 0, and consider the nonlinear parabolic partial differential equation 
uz — (D(u)uz), = ®(2,t), (x,t) € 2 = (0,1) x (0,7), (6) 
with initial condition 
u(x,0) = s(a), x € [0,1], (7) 
and boundary conditions 


u(0,t) = f(t), t € [0,7], (8) 
u(1,t) = A(t), t€ [0,7], (9) 


where D(u) = a(t)u+ W(t) > 0, anda, b, s and f are known functions 
such that, b(t) is far from zero in (0, T]. 

If h is given, then the problem (6)-(9) is a direct problem which is 
solvable by means of common numerical and approximation approaches for 
solving PDEs, such as finite difference method [7], finite element method 
[7], radial basis functions [4], homotopy perturbation method [13], Adomian 
decomposition method [1] and so on. 

Now, suppose that h is unknown. Then the problem (6)-(9) becomes an 
inverse problem. Consequently, an overspecified condition, such as 


uz(0,t) = g(t), t € [0,7], (10) 


where g is a known function, is used. 

For positive integer n, let At = i = oan t; =jAt,j € Irn = {1,2,..., n}. 
Put ug(a) = u(x, 0) = s(x), a; = a(t;), 6; = d(t,;), ©;(@) = O(2,£,), for 
any j € Jz, such that they are given fixed nodes. Similarly, we consider the 
u,j(z), as the approximated value of u(z,t;), 7 € Jn. 

Using the backward finite difference scheme for the term u; in the form 


u(x) ~ k(u;(x) — uz-1(2) ), jEdIn, 


and substituting in Eq. (6), a system of second order ordinary differential 
equations with respect to x is obtained. We have 


B(uj(2) ~ uya(2)) ~ Ff (alts ug(2) + B(ts)) Sry @)} = Oj(@), <7 Som, 


or 
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=] 

—~ ©, (x). 11 

b(t;) 5 (2) ( ) 
For simplicity, define u(x) = (ui(z), ue(z),...,Un(x))?. Then we can 


write Eq. (6) as follows 
Au = L,u— Nu= WV(z,t), 


where W(x, t) = (age. Ee) Moreover, L, and N are the linear 


and nonlinear parts of the operator A, respectively, and are as follows: 


d2 
co dues 
d 
Nu= -Ma 7 (D(u,u)) +Myu+m, 
v 
where ; 
b ti) 0 
Dis) Bia) 
oa 2 2 
m= Cas s(x), 0, 0)" Mi =k i ’ 
0 TOs) Bia) 
; a(t a(tn 
Me = diag ( b(t)? nae ’ ea i 
and 
D(u(x), v(2)) = (ui (x) v1 (x), uo(x) Fe v2(a), aes Un (2) Fn (x))”. 


After twice integration of Eq. (6) with respect to x, and applying the condi- 
tions (7)- (9), we obtain 


u(a)-og-t- ff Nu (2) ards = [ [ weaeae 


where g=(g(ti), --- , g(tn) ), and f=( f(t), --. , f(tn) )*. 
Now, using HPM and [8,26], we choose a convex homotopy such that 


H (v(x), p) = v(x) — h(a) — of A Nv(«) dxdx = 0, (12) 


and 


where 
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h(a) )=eg+t+ [ [ ) dxdax, 
and v(x) =((2,t1), ... v(a,tn) ) . Furthermore, Eq. (12) gives 


v(x) = h(z) +of a Nv(a) dadz. (13) 


Combining Eq.s (11) and (16), we obtain the following results 


v(x )=ag+f+fp ya ) dada + 


EL on tw ere rn (14) 


where u,(x) = (s(7),uy(a),...,Un—1(x))7. Thus, it is concluded that 
vo(x) = = ) =eettt [ [ow x)dadx = (vo(2,t1),--+ ,vo(2,tn))” , (15) 
and 


= [f° {xn cote) - ao) ~ Ma £ D(vo(a),vo(a)) dade 
eer ye (16) 


= (vi (2, t1), 
The above relations are obtained by equating the terms with identical powers 
of p in Eq. (19). The approximate solution is 


u(x) © vo(xz) + vi(x) = v(x). (17) 


4 Convergence and stability analysis 


In this section, we use continuity of u(a,t) on the compact domain Q), 
and prove that v(x) ~ vo(z) + vi(x) depends continuously on the data. 
Therefore, adding a perturbation term to a, b, f, g and ®, an upper bound 
for errors of their solutions are found. In each case, we show that as the 
perturbation term tends to zero, the solutions errors vanish. 


Lemma 1. Let M(t =fy |®(x,t)|dx > 0 is a bounded function such that 
M(t) <M for any : <t<T, and w(x) correspond to v(x), where b(t) 
is perturbed by db(t) in Equation (6). Then we have 


|dd(t;)| 
[b(E,) (Os) + SBC)” 


O89) “HS [0el, 


|Go(x, tj) — vo(z, ty)| = 
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consequently, if |db(t)| > 0, then |do(a,t) — vo(a,t)| > 0, for any t = 


tj, j=l, +++, n, and O<a<l. 


Proof. Using the Equation (15), we have 


vo(x,t;) = f(t;) + g(t mol [a «dade. (18) 


If b(t;) is perturbed by 6b(t;), then 


folety) = It) +204) — aE a. [a n)dede. (19) 


Now from Equations (18) and (19), we obtain 


vo(a, t;) = 00 (x, tj) 5 b(t; ) (b(t or va i ie [a nee 


consequently 


[db(t5)| 
[b(t )(b(t;) + db(t;)) 


and this completes the proof. 


|vo(a, t3) — So(x,t;)| < M, foranyO<j<n, 


Lemma 2. Let M(t) and M_ are as defined in lemma 1, and %1(z,t;) 
is the value of vi(x,t;) for anyt=t;, jg =1, +: , n, when d(t;) is 
perturbed by b(t;) + S(t j) as in problem Equation (6), such that 6b(to) = 0. 
If |db(t;)| + 0, then ens — vi(a,t;)| 3 0, for any O< «<1, and 
l<jg<n. 


Proof. Similar to detailed proof presented for lemma 1, assume 


(aes -[ Ge CHCS OC) 


 {vo(2, ts) vol, t;)})}dadz. 


Suppose that b(t;) is replaced by Ob(t,;) + db(t;), for any 7 = 1, --- , n. 
Then we have 


(2,45) =f fe ae 7 


a(t;) d 
b(t;) + o() ‘der aon tj 5) Gy 00 o(ax,t;)})}dada, 


or 
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i}. 1 A es 

r \Go(x,t;) — vo(a,t;)| | |aj-1(@) — uj—r(z)| 

01(a,t;) — vi(a,t;)| <k ' 

pea UAE) To) + o0(E) [b(t;) + 5b) 


|db(t;)|C oa 
DE )(6(E;) + 3B I 

la(t;)||db(t;)| M 
b(t; )(D(t;) + 6b(t;))| [b(t )(b(t;) + 6b(E,))| 


x max {|d0(z, t3)|, |vo(x, t3)|} 
+ IF (ts) 9(ts)| 
t;)| 


when C' € R is an upper bound for |vo(,t;) — uj;—1(x)|. Therefore, for fixed 
k, from lemma 1, it is derived, 


gHlye yon. 
_ leo ti) + FG 
2|( 


), 


li v t;) — t;)| —>0 
Pes tae tac sh ties) | : 


for O<a<1, and l<j<n. 


Remark 1. By an induction, it is shown that, when |db(t;)| —> 0, the 
second term of the integral in the above inequality vanishes. 


Theorem 1. Suppose that v(zx,t;), vo(z,t;), do(a,t;), vi(x,t;), O1(%,t;), 
M, M(t), b(t;) and db(t;), are the same as defined in lemmas 1 and 2. Then 
u(x, t;) = vo(x,tj) + vi(x,t;) depends continuously on the data. 


Proof. Obviously, by considering lemmas 1 and 2, the statement of Theorem 
1 is proved. 


Theorem 2. Suppose that da(t;) is the perturbation term that perturbs a(t;) 
to a(t;) + da(t;), and vo, vi, to, t1, M, M(t), b(t;) and db(t;) satisfy 
assumptions of Theorem 1. Then v(x,t;) depends continuously on the data. 
Proof. The first part of u(x,t) is independent of a(t). Then, using Theorem 


1, there is nothing to prove for vo(x,t;), 1 <j <n. Now we just prove that 
v1(z,t;) depends continuously on the data. We have 


A _ [” [%,, bola, t3) — Gj-1(@) ~~ a(tj) + da(t;) , d 
ne) = f 5 UE) EEG) BE) SBME) ae 


{60 (2, t)) -t0(0,t;)}) bande, 


then, it is concluded that 


jatet)- ates fo f Pei eel 4 Bee 
(3) 


tj) + db(t;)| |b(t;) + db(t,)| 
\dadx 


[b(t )(b(t;) + 68(¢;))| 
a(t; )|]66(¢5)| M 


[b(t )(b(ty) + 6b(E;))| “To(E YUE) + db(t;))| 
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nas {ole t)h lone t0 


vale ty) + Liles)? | + |F(ts)9(ts)| 
2|b(t;)| 
L Jsa(t,)| 20etT 


) 


Similarly, it is shown that, the approximate solution v(x) in (17), depends 
continuously on the data, when f(t), g(t) and ®(z,t) are perturbed by small 
perturbation terms in their domains. So, we give the following theorem. 


Theorem 3. Let f(t), g(t) and ®(x,t) be the same as defined in Equation 
(6). If 
f(t) > ft) + 6f(), 
g(t) —> g(t) + dg(t), 
®(x,t) +> ®(2,t) + 6®(z,t), 


and v'+— v + dv, then |du(a,t)| —> 0, when 


max {|6f(¢)| + |dg(¢)| + |6®(2, t)|} — 0. 


O<t<T 


Furthermore, we have 


ele.) — vlet)1 < [bre + ldate] + Se 
xj "| i yh ei {|0(x, tj) — vo(a, t;)| + |aj-1 — uj—i|}dadx 


2\b (t; rami 06(a, tj) v6 (a, ty)]- 


Proof. We have 


Holaty) = H(t) + OFC) + 20lts) + 2806) — Ges ff Corts) + 40(e,1)))drde, 
J 
and so 


P(x, ts) 


Go (a, t;) — vo(a,ts)| < \6F(ts)| + 6g(ty)| + 2o(t;)| 


(20) 


So, |éo(x, t) — vo(a,t)| —> 0, when 


max {df (¢)| + |6g(¢)| + |6®(2, t)|} — 0. 


O<t<T 
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Similarly, putting 

ete ee al ree he, CALE wy ee cee ee . 
artestyy= ff {ae Ceolests) ~ 0)—1(0)) ~ DF (oo(e, ty) Z 02, 45))} dd 


which via lemma 2 and Equation (20), simplifies to the form 


: k cant aaeee , 
mee ie) Se ff tieote.ts) = voles ts)| + [aya ~ ujal} dee 
ine o Jo 
a(t;)| \ 2 2 
t;)- t;)|. 21 
BC Pts) ~ BC) (21) 


Now, |é1(x, t) — v1(x, t)| —> 0, whenever 


max {|df(t)| + |d9@)| + [S®(@, t)|} — 0. 


By adding the two sides of Equations (20) and (21), we obtain 


i |5®(x, t,)II 
|o(x, t3) — v(a, t3)| ped + |dg(t;)| + “bE 
ca in | ‘i {|60(x, t;) — vo(x, t3)| + |aj—-1 — uj—i|}dadx 


la(ts)| ) 02 2 
a 2|b(t;)| |0o (x, tj) 35 vo (x, t;)|- 


Finally, |(a, t) — v(x, t)| —> 0, whenever 


max {df (¢)| + |6g(¢)| + |6®(2, t)|} — 0. 


O<t<T 


In the next section, a necessary condition for convergence of the approx- 
imate solution, when the step size At, tends to zero is obtainsd. 


5 Convergence conditions for the problem (6)-(9) 


In this section, we apply the error term of finite difference method in the 
relations (15), (16) and (17) and then convergence condition of the solution 
will be investigated. So, first, we give the following theorem. 


Theorem 4. Let |At b(t)| > 1 for anyO < t < T. If |At]| —> 0, then 
uj (v7) —> u(az,t), for anyO<j<n, O<a<1. 


Proof. Using Taylor’s series expansion, we have 
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u;(a) ~ u(x,t) + 6®, 


where 60 = At24#@9) and tj < 0; < ty41. 
So, from (18) we obtain 


vo(z,t;) = f(t;) + xg(t; mol [ a x)dxdx, 
or 


v0j = f (tj) + vg(t; mo, [ast x)dada + 5®;. 


Clearly, if 66; —>0, then vo; —> vo(z,t;). 
Again, by (2 ) we have 


a= a BE) 0 x,tj) — uj—1(#)) — MH) (4 tyo(2,t)) vole, t;)})dada. 


b(t;) 
thus 
x zx k 
v= 7 b(t) i ep tley SOD a4) 
b(E,) ‘aa 6% bj) 7 (v0j — 5b;)} dade 
+ 68;, 
and 
O15 = 01 (2,1;) one ie (d@; — 6®,_,)dadx 
alts) ae; i 6®.d 6® d. 
Xi) S| 04 j c- [ 5-05 x}. 


That means if 66; —> 0, then v,; —> vi (2,t;), and v; — v(2,t;). 


6 Numerical results 


In this section, we give a numerical example. 
Let 


Ou, 7 


{(= e ‘u+ (t+5) ee = age (x,t) € [0,1] x [0, 1], 


u(z,0)=27, 0<a#<1, 
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u(0, t) = t, <t<l, 
u,(0,t)=0, O<t<1. (22) 
Obviously 
®(x,t) = —it-9, =. b(t) = (+5) e*. 
The exact solution is u(x,t) = x?e' +t. We obtain the approximate solution 
by applying equations (15), (16) and (17), at = 0.1, 0.2,... , 1, where t = 


0.25, 0.5, 0.75, 1, and we assume that At = 0.25. Consequently the solution 
will be constructed in the form 


2 
Vo; (x) = Aw, t;) = tj ea ( tt 9) - ’ 


wes { Aay(vo, (2) Ge te 4)) = eae ( vo, (x) vo, (zx)) dxdz, 


for 7 = 1,2,3,4. 

The exact solution, approximate solution and relative error for the above 
problem are given in Tables 2—5 and6—9att=t,;=jAt, j =1,2,3,4. 
To illustrate stability, according to Table 1, we enter some noise terms into 
data functions in Eq. (22). 


TOIT . 
1.0 nye AE + 


Figure 1: Approximate (---) and exact solution of u(z,t;) in t = 0.25 and t = 0.5. 
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Figure 2: 
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Table 1: Perturbation terms in problem (22) 


-) and exact solution of u(x,t;) int = 0.75 and t = 1. 
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a5 


bs 


fi 


I5 


®; 


FP wN Blo, 


-0.006486161986 
0.01848017099 
-0.005506853028 
-0.006487155975 


-0.001002685512 

-0.0008600955762 
0.002919491298 
-0.001056710212 


0.001219370604 
-0.0002264198307 
-0.0004964518833 
-0.0004964988898 


-0.001002685512 
-0.0008600955762 
0.002919491298 
-0.001056710212 


0.001219370604 
-0.0002264198307 
-0.0004964518833 
-0.0004964988898 


Table 2: Exact and approximate solution of uj(x) at t; = 0.25 with and without pertur- 


bation terms in a and b 


exact solution 


relative error 


perturbed solution 


relative error 


x approximate solution 
0.1 0.2628402542 0.2628483725 3.08 x 107° 0.2628562616 6.09 x 10-5 
0.2 0.3013610167 0.3013841491 7.67 x 107° 0.3014170267 1.85 x 10-4 
0.3 0.3655622875 0.3655793088 4.65 x 107° 0.3656582382 2.62 x 1074 
0.4 0.4554440667 0.4553871491 1.24 x 10-4 0.4555397988 2.10 x 10-4 
0.5 0.5710063542 0.5707422861 4.62 x 10-4 0.5710055740 1.36 x 10-6 
0.6 0.7122491501 0.7115606556 9.66 x 10-4 0.7119833904 3.37 x 1074 
0.7 0.8791724543 0.8777395127 1.62 x 1073 0.8783850345 8.95 x 1074 
0.8 1.071776267 1.069157431 2.44 x 107% 1.070106259 1.55 x 1073 
0.9 1.290060588 1.285674303 3.40 x 10-3 1.287026770 2.35 x 10-3 

1 1.534025417 1.527131339 4.49 x 1073 1.529010242 3.26 x 1073 
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Table 3: Exact and approximate solution of uj;(x) at t; = 0.5 with and without pertur- 


bation terms in a and b 


exact solution 


approximate solution 


relative error 


perturbed solution 


relative error 


0.5164872127 
0.5659488508 
0.6483849144 
0.7637954034 
0.9121803178 
1.093539658 
1.307873423 
1.555181613 
1.835464230 
2.148721271 


0.5165050073 
0.5660064728 
0.6484638498 
0.7638099282 
0.9119513987 
1.092769640 
1.306121720 
1.551841629 
1.829741743 
2.139614496 


3.44 x 10-5 
1.01 x 10-4 
1.21 x 10-4 
1.90 x 107° 
2.50 x 10-4 
7.04 x 10-4 
1.33 x 1073 
1.09 x 10-3 
3.11 x 10-3 
4.23 x 1073 


0.5164663253 
0.5658439475 
0.6480689075 
0.7630349271 
0.9105938906 
1.090556514 
1.302693219 
1.546735182 
1.822375591 
2.129271087 


4.04 x 1075 
1.85 x 10-4 
4.87 x 10-4 
9.95 x 10-4 
1.73 x 1073 
2.72 x 107 
3.96 x 1073 
5.43 x 1073 
7.13 x 107 
9.05 x 10-3 


Table 4: Exact and approximate solution of uj;(x) at tj = 0.75 with and without pertur- 


bation terms in a and b 


exact solution 


approximate solution 


relative error 


perturbed solution 


relative error 


0.7711700002 

0.8346800007 

0.9405300015 
1.088720003 
1.2779250004 
1.512120006 
1.787330008 
2.104880011 
2.464770014 
2.867000017 


0.7711983689 

0.8347737522 

0.9406671824 
1.088781041 
1.278980053 
1.511092655 
1.784912741 
2.100201772 
2.456691211 
2.854085288 


3.36 x 107° 
1.12 x 10-4 
1.45 x 10-4 
5.60 x 10-5 
2.11 x 10-4 
6.79 x 10-4 
1.35 x 1073 
2.22 x 10-3 
3.27 x 10-3 
4.50 x 1073 


0.7712065247 

0.8348115589 

0.9407717933 
1.089016095 
1.279447063 
1.511943216 
1.786361766 
2.102541999 
2.460309231 
2.859479316 


4.73 x 10-5 
1.57 x 10-4 
2.57 x 10-4 
2.71 x 1074 
1.54 x 10-4 
1.16 x 10-4 
5.41 x 10-4 
1.11 x 10-3 
1.80 x 10-3 
2.62 x 10-3 
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Table 5: Exact and approximate solution of u;(x) at tj = 1 with and without perturba- 


tion terms in a and b 


xX 


exact solution 


approximate solution 


relative error 


perturbed solution 


relative error 


1.027182818 
1.108731273 
1.244645364 
1.434925092 
1.679570457 
1.978581458 
2.331958096 
2.739700370 
3.201808281 
3.718281828 


1.027222369 
1.108860738 
1.244829262 
1.434986046 
1.679134916 
1.977027872 
2.328368191 
2.732814161 
3.189983398 
3.699457703 


3.85 x 1075 
1.16 x 10-4 
1.47 x 10-4 
2.42 x 10-5 
2.59 x 10-4 
7.85 x 10-4 
1.53 x 10-3 
2.51 x 10-3 
3.69 x 1073 
5.06 x 10-3 


1.027311557 
1.109229081 
1.245701389 
1.436644053 
1.681940678 
1.981444677 
2.334981455 
2.74235 1004 
3.203330804 
3.717678990 


1.25 x 10-4 
4.48 x 10-4 
8.48 x 10-4 
1.19 x 10-3 
1.41 x 1073 
1.44 x 1073 
1.29 x 1073 
9.67 x 10-4 
4.75 x 10-4 
1.62 x 10-4 


Table 6: Exact and approximate solution of uj;(x) at t; = 0.25 with and without pertur- 


bation terms in f, g and ® 


exact solution 


approximate solution 


relative error 


perturbed solution 


relative error 


0.8 
0.9 
1 


0.2628402542 
0.3013610167 
0.3655622875 
0.4554440667 
0.5710063542 
0.7122491501 
0.8791724543 
1.071776267 
1.290060588 
1.534025417 


0.2628483725 
0.3013841491 
0.3655793088 
0.4553871491 
0.5707422861 
0.7115606556 
0.8777395127 
1.069157431 
1.285674303 
1.527131339 


3.08 x 107° 
7.67 x 107° 
4.65 x 107° 
1.24 x 10-4 
4.62 x 10-4 
9.66 x 10-4 
1.62 x 10-3 
2.44 x 1073 
3.40 x 10-3 
4.49 x 1078 


0.2630645568 
0.3015988237 
0.3657914404 
0.4555956649 
0.5709460580 
0.7117584843 
0.8779301112 
1.069339409 
1.285846152 
1.527291415 


8.53 x 10-4 
7.89 x 10-4 
6.26 x 10-4 
3.32 x 1074 
1.05 x 10-4 
6.88 x 10-4 
1.41 x 1078 
2.27 x 10-3 
3.26 x 1073 
4.38 x 1073 
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Table 7: Exact and approximate solution of uj;(x) at t; = 0.5 with and without pertur- 


bation terms in f, g and ® 


exact solution 


approximate solution 


relative error 


perturbed solution 


relative error 


0.5164872127 
0.5659488508 
0.6483849144 
0.7637954034 
0.9121803178 
1.093539658 
1.307873423 
1.555181613 
1.835464230 
2.148721271 


0.5165050073 
0.5660064728 
0.6484638498 
0.7638099282 
0.9119513987 
1.092769640 
1.306121720 
1.551841629 
1.829741743 
2.139614496 


3.44 x 10-5 
1.01 x 10-4 
1.21 x 10-4 
1.90 x 107° 
2.50 x 10-4 
7.04 x 10-4 
1.33 x 1073 
1.09 x 10-3 
3.11 x 10-3 
4.23 x 1073 


0.5154115152 
0.5648920604 
0.6473145948 
0.7626119527 
0.9106908844 
1.091432842 
1.304694990 
1.550311436 
1.828094690 
2.137837344 


2.08 x 1078 
1.86 x 1073 
1.65 x 1073 
1.54 x 1073 
1.63 x 10-3 
1.92 x 10-3 
1.30 x 1073 
3.13 x 1073 
4.01 x 1073 
5.06 x 1073 


Table 8: Exact and approximate solution of uj;(x) at tj = 0.75 with and without pertur- 


bation terms in f, g and ® 


exact solution 


approximate solution 


relative error 


perturbed solution 


relative error 


0.7711700002 

0.8346800007 

0.9405300015 
1.088720003 
1.2779250004 
1.512120006 
1.787330008 
2.104880011 
2.464770014 
2.867000017 


0.7711983689 

0.8347737522 

0.9406671824 
1.088781041 
1.278980053 
1.511092655 
1.784912741 
2.100201772 
2.456691211 
2.854085288 


3.36 x 107° 
1.12 x 10-4 
1.45 x 10-4 
5.60 x 10-5 
2.11 x 10-4 
6.79 x 10-4 
1.35 x 1073 
2.22 x 10-3 
3.27 x 10-3 
4.50 x 1073 


0.7736467637 

0.8372983250 

0.9433190828 
1.091611970 
1.282042479 
1.514440029 
1.788599719 
2.104284426 
2.461227246 
2.859134249 


3.21 x 10-3 
3.13 x 1073 
2.96 x 10-3 
2.65 x 10-3 
2.18 x 10-3 
1.53 x 1073 
7.10 x 10-4 
2.82 x 10-4 
1.43 x 1073 
2.74 x 107 
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Table 9: Exact and approximate solution of u;(x) at tj = 1 with and without perturba- 


tion terms in f, g and ® 


x exact solution 


approximate solution 


relative error 


perturbed solution 


relative error 


0.1 1.027182818 1.027222369 3.85 x 107° 1.025635295 1.50 x 
0.2 1.108731273 1.108860738 1.16 x 10-4 1.107171622 1.40 x 
0.3 1.244645364 1.244829262 1.47 x 10-4 1.242968576 1.34 x 
0.4 1.434925092 1.434986046 2.42 x 107° 1.432882007 1.42 x 
0.5 1.679570457 1.679134916 2.59 x 1074 1.676712572 1.70 x 
0.6 1.978581458 1.977027872 7.85 x 1074 1.974208173 2.21 x 
0.7 2.331958096 2.328368191 1.53 x 10-3 2.325067038 2.95 x 
0.8 2.739700370 2.732814161 2.51 x 107% 2.728941436 3.92 x 
0.9 3.201808281 3.189983398 3.69 x 107% 3.185441974 5.11 x 

1 3.718281828 3.699457703 5.06 x 107% 3.694142422 6.49 x 


10-3 
10-3 
10-3 
10-3 
10-3 
10-3 
10-3 
10-3 
10-3 
10-3 


7 Conclusions 


The HPM for the one-dimensional inverse problems has been presented. The 
method described is mathematically simple and computationally effective. 
As we see in Tables 2 — 9, small errors in the data make small errors in 
the solution, so that, the solution depends continuously on the data. In 
this paper, the noise terms that are shown in Table 1, are made randomized 
and have standard normal distributions. Maple 16 packages have been used 
to compute the solution before and after adding noise terms. Rapidity , 
accuracy and stability are advantages of this formulation. 
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